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Abstract 

In this contribution we apply the QCD dipole picture combined with A^-factorization 
to get predictions for deep-inelastic scattering on an onium. Assuming renormalization- 
group factorization, we get predictions for the F2,Fq and R = Fl/Ft proton structure 
functions. We obtain a three-parameter fit of the 1994 HI data in the low-x^-, moderate- 
Q 2 range. Fg/F2 and R are then predicted without any additionnal parameter. The 
BFKL dynamics contained in the dipole model is shown to provide a relevant model in 
describing the HERA data. The prediction for F2 and Fq are compatible with next-to- 
leading order DGLAP analysis. By contrast, R is expected to be much lower at small 

Xbj- 
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1 Introduction 



In this contribution, we show that the QCD dipole picture [|1], 0, |3|, f|], which contains the BFKL 
dynamics, is a pertinent model for describing the proton structure functions at HERA in the 
\ow-Xbj and moderate-Q 2 range. The most recent HERA 1994 published data || provide a 
possibility of distinguishing between the different QCD models for small-Xfej physics. Besides 
the DGLAP evolution equations || based on the renormalization group evolution, the BFKL 
hard Pomeron |7j should be relevant at small Xbj- It is thus a challenge to get observables which 
could distinguish between these two different theoretical predictions, and to compare these 
observables to the data. The perturbative dipole model, combined with the /^-factorization 
(or Regge factorization) ||, allows us to implement the BFKL dynamics into deep inelastic 
scattering. We then get predictions for F 2 , Fq and R = Fl/Ft proton structure functions. 



2 BFKL dynamics in the QCD dipole framework 

In this section, we perform an analyse of e ± — p deep inelastic scattering at low x^, based on 
fc T -factorization and dipole color model, as illustrated in figure [p. 




Figure 1: /c^-factorization and dipole model applied to e — p deep inelastic scattering. 



In the Regge limit, one can apply the /^-factorization tool pi in order to extract a photon of 
virtuality Q 2 off a proton. It involves the elementary Born cross-section a ig /Q 2 of the process 
7 g{k) — > q q. Here the gluon is off-shell, quasi transverse, with a virtuality k 2 ~ A; 2 . One also 
has to introduce the unintegrated gluon distribution density at a factorization scale Ql, which 



-I 



is related to the usual gluon distribution by 

G(x bj , Q 2 , Ql) = [ Q d 2 kf(x bj ,k,Ql). (2.1) 

<j 

We first deal with a dipole of transverse size Xoi, which can be part of an heavy onium (i.e. 
an heavy qq pair) or extracted from a proton as will be emphazied later. The fc^-factorization 
reads, for the total 7* — dipole cross-section, 

Q 2 tf*(x bj ,Q 2 ;x 2 01 ) = J d 2 kj^ ^ a^ g (x bj /z, T(z, k; Xqi). (2-2) 

But we have now to couple the gluon to the softest dipole which arise in the cascade. This can 
be fulfilled by using a second fc-r-factorization. It involves the elementary Born cross-section 
o^jk 2 of the process d(x) g{k) — > d(x) for a dipole of transverse size x and a soft gluon of 
virtuality fc 2 . This /cy-factorization can be expressed by 

i 2 ^"(^,i;a;oi) = / — P / ^r n (x 01 ,x,z')d- ld (z/z',x 2 k 2 ), (2.3) 
J xf Jo z' 

where n(x 01 ,x, z) is the distribution density of dipoles of transverse size x with the smallest 
light-cone momentum in the pair equal to zp + in a dipole of transverse of total 

momentum p + . Computing <r 7d with an eikonal coupling of the gluon ||, one gets 

Q 2 c*.(x bj ,Q-) = J d 2 kJ Q ^-a ig (x/z,=p) J n(xQ U x, z) 

x4vra ^ (2 - - e~^)\. (2.4) 
(2vr) 4 v J k 2 K ' 

Following [|IJ, we now introduce the double Mellin-transform of n with respect to the rapidity 
Y = \wz\jz and the transverse size x {z\p + is the light-cone momentum of the quark part of 
the dipole (x 01 )): 

n(x 01 ,x,Y) = J ^e wy r^(x 01 ,x). (2.5) 

and 2 

n^xjox^x) = J ^ (^) n w (7) (2.6) 

Note that (|2.6|) is a decomposition into conformal three-points correlation functions, due to the 
conformal invariance of the dipole integral kernel. One then gets |l| 

"-(7) = Jc F ( s ( 2 - 7 ) 

where 

x ( 7 ) = ¥(l)-I*( 7 )_Itt(l- 7 ) and *( 7 ) = ^£. (2.8) 

o ig has been calculated for different polarizations of the incoming photon in ||. We introduce 
the corresponding double Mellin transform of this cross-section, namely 



(2.9) 



or equivalently 

l 2 \ 2 f d 1 ( I 2 V 7 1 



= m ' ^\q~ 2 7 M7) - (2 - 10) 



The expression for a^, now reads 

g a r K-,Q ;x M ) = 4vr / d k J — J — J —e (_j 



The integration with respect to the polar angle of x can easily be performed and leads to a 
Bessel function. One has then to integrate over x, namely 

/ _ 4vr (l _ J (kx)) — = 4vr(A;xoi) 27 ■ 7 = 47r(^ i)%( 7 ). (2.12) 

•/ x Vxoi/ 7 1(1 + 7) 

The integration over /c gives 7 = 7'. Since n vw (formula (|2.7| )) exhibits a pole at uo p = 9 ^ £ -x{l)i 
the uj integral finally yields 

Q2 4(^,Q^ Sl ) = ^/^ ft „,( 7 )^Ma iQ ^ e ^*)^. (2.13) 



47r 2 a e . m rv ' Ui/ 1 i 2m " pv,/ 7 

In the following we neglect the dependence of h Up {^) on u p ||. The initial dipole state is 
supposed to be well localized in transverse space, and its transverse size gives a perturbative 
scale. We then get for the total cross-section 7* — dipole(x m ) and for the related structure 
function: 

Q 2 „d „2. ji n _ rtf- ^2. j2 ^ _ 2 ® N c f d 7 ,^ 2 _ a ^^(7)^^)1^ 



-a d r (x bj} Q 2 ;x 2 01 )=F*(x bj ,Q 2 ;x 2 01 ) = f^l [ |1 (Q^^t)- 



4vr 2 a em 7 v 7 uu TV ' U1 ' 7T 7 2m ^ Ui/ w ' 7 

(2.14) 

In order to deal with deep inelastic scattering on a proton, we will suppose that inside the 
proton there exist configurations of dipole type revealed by the high-energy process (uu for 
example) which are well transversally localised, with a typical transverse scale Xq X . The dis- 
tribution w(j, Xqi, Qq) of such configurations cannot be computed perturbatively, but, using 
renormalisation group properties, this quantity is expected to have the following behaviour (see 



Toll for a similar analysis applied to the leading order QCD evolution equation): 

w( 1 ,x 2 01 ;Ql)=w( 1 )(xl 1 Ql)-\ (2.15) 

Ql is a scale which is typical of the proton and expected to be non-perturbative. We get the 
final result: 



— X(7)ln— V (2.16) 



At this stage some comments are in order: i)We insist on the point that these initial pertubative 
configurations, even if they cannot be described easily, should in principle be present inside the 
proton, and their distribution non-negligible (in contrary to cc pairs for example, which, from 
J/ty production, are known to be less relevant). Their evolution over a large range of rapidity 
is responsible for the cascade which contains the physics of the BFKL Pomeron. ii)This BFKL 
dynamics can also be implemented when using the original BFKL kernel combined with Regge 



factorization. It involves an initial unintegrated gluonic distribution, which is then evolved 
to lower values of x^ by the cascade of reggeized gluons. Here again, this initial distribution 
cannot be predicted and has to be put by hand [OH. 

We now apply formula (|2.16|) to peculiar structure functions, namely 
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(2.17) 



where -FV(l) is the structure function corresponding to transverse (longitudinal) photons and 
Fq the gluon structure function. The coefficient functions 



h T \ a (r(i- 7 )r(i + 7 )) 3 1 / (1 + 7)(1- 

h L J 3tt 7 T(2 - 2 7 )T(2 + 2 7 ) 1 - § 7 \ 7(1-7) 



(2.18) 



were computed in ref [Q . The coupling constant a enters the coupling of the virtual photon to 
the off-shell gluon, and thus in formula CP .18 ). It also enters the effective coupling of the dipole 
kernel. We shall take it identical in both cases. We also keep a non running, even if it might be 
possible to implement the running of at in the dipole cascade [ll|. Indeed, this would involve 
next-to-leading terms which are beyond this model and not already known. Nethertheless, the 
running of a is not expected to drastically change the shape of the dipole cascade. 
When is small, the 7— integration can be performed by the steepest-descent method. Consid- 



ering w(j) as smooth and regular near 7 
anomalous dimension 7 - 



2' 



we expand the x function around the asymptotic 



, . We obtain a saddle point at 



7, 
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where 
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The approximation of expanding xil) around \ is valid when 
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that is the small Xbj, moderate Q/Qo kinematical domain. This yields 



-1 



(2.19) 



(2.20) 



F2 = Ft + Fl 



r< 1/2 Q 

Ca ' — exp 
Qo 



(2.21) 



where 



«p — 1 



4aN c In 2 



7T 



(2.22) 



Thus, F 2 depends only on 3 parameters, C, Qo and otp. Fitting F 2 with this form, one can get 
a prediction for F G and R = F L /F T . Namely, 



Fg 
F 2 



h T + h L 
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r(2 - 2 7s )r(2 + 2 7s 



and 
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H 2 (r(i-7 s )r(i + 7s)) 3 



—h 



ls{l -Is) 



(2.23) 



(2.24) 



(1 + 7 S )(1-1)' 

where 7 S is given by the expression ( |2.19| ). Note that the overall non-perturbative normalization 
C does not enter R and Fq/F 2 . 



1 



3 F 2 fit and prediction for Fq and R 



A fit of the 1994 HI data || has been performed. We have kept only bins in Q 2 with Q 2 < 150 
GeV 2 , in order to stay in the range where the QCD dipole model is relevant and to insure that 
our saddle point method can be trusted (see the limit ( |2.2Up ). This point is discussed in ref. 
fL9f . Note that our fit does not depend strongly on that cut value. The fit parameters are C, 
Qo and ap. The obtained fit for F 2 is displayed in figure || (the highest Q 2 point at 5000 GeV 2 
is not displayed for convenience). We get a \ 2 °f 101 for 130 points, and the following fitted 
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Figure 2: Comparison of the fit of the HI data for Q 2 < 150 GeV 2 with all the 1994 HI data. 
The discrepancy at high Xbj and high Q 2 is outside the validity of our model. 

parameters: 

a P = 1.282 

Qo = 0-627 GeV 



C = 0.077. 



(3.25) 



The obtained value for the hard Pomeron intercept ap is in agreement with other determinations 
applying BFKL dynamics to F2 data at HERA [[11], |5|. The corresponding effective coupling 
constant is a ~ .11, close to a(M Zo ) used in the HI QCD study. However, the obtained value 
for the Pomeron intercept is relatively smaller than the value one would expect («p ~ 1.5), and 
the corresponding effective coupling is surprisingly small. This might be due to next-to-leading 
order terms, responsible for this effective decrease of a p. This is indeed the case for example 
when one attempt to sum up both soft and collinear singularities and if one takes into account 
energy-momentum conservation [TB|, [TJ|. The value of Qo corresponds to a transverse size of 
.3 fm, which is, as expected, a non-perturbative scale typical of deep inelastic scattering on a 
proton target. The non predictable parameter C is linked to the non-perturbative input w(^). 
It is quite remarkable that one can well reproduce the data in the low-Q 2 range with only 3 
parameters. It confirms a previous study [[3j| using HI and ZEUS 1993 data | 16| . Deviations 
from the fit at higher-Q 2 values have several origins. First, the valence contribution is not 




Figure 3: Prediction for the gluon density F G (xbj, Q 2 ) (full line) compared with the HI predic- 
tion based on their next-to-leading DGLAP QCD fit. In dotted line the one loop approximation 
is displayed. 

contained in the dipole model, and it is well known that it should dominate at higher-x^- 



values. Because of the HERA kinematical constraints, this corresponds to bins at higher-Q 2 
values. Second, we used an expansion of the dipole eigenvalue around ~, which is valid only for 
moderate-Q 2 values. A phenomenological input for the valence contribution and a numerical 
treatment of the 7 integral is expected to improve the fit of the HERA data in the higher Q 2 
range. 

Relation ( |2.23| ) provides now a precise determination of the gluon density. We use for Q 
and a the values obtained from the F 2 fit. In figure |3] we show the comparison between the 
published HI determination of the gluon density (grey band), based on a next-to- leading order 
DGLAP analysis, and our prediction (full line). We also display the prediction obtained from 
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Figure 4: Prediction for the R (full line) compared to the Altarelli-Martinelli prediction used 
in the HI study. In dotted line the one loop approximation is displayed. 

first order perturbative expansion of the coefficient functions h, in order to see the resummation 
effects. We notice that for x < 10~ 2 , our prediction is in the range quoted by HI collaboration 
||. Thus, the gluon distribution does not provide a good tool in distinguishing between DGLAP 



and BFKL dynamics. 

The prediction for R is rather different. In figure |], we display the Altarelli-Martinelli 
l7j prediction for R, used in the HI study f|, in comparison with our prediction. We also 



show up the resummation effect for h. The discrepancy between AM and BFKL predictions is 
rather strong. In particular, the limiting value of the dipole model is 2/9, that is much lower 
than the AM prediction [M . R could thus provide a good observable in distinguishing between 
the two possible evolutions. 



4 Conclusion 



In this contribution we have shown that combining the dipole model with fcr-factorization, it is 
possible to implement the BFKL dynamic in order to describe the proton structure functions. 
With only 3 parameters, the most recent 1994 HI data can be well described in the range 
Q 2 < 150 GeV 2 . We obtain a prediction for Fq which is close to the one based on next-to- 
leading order DGLAP equations. In contrast, the prediction for R is significantly different. 
Measurement of this quantity would thus be of great interest in the low-x^- range. 
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